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1 Introduction 

A wide range of phenomena observed in heavy-ion collisions at RHIC and LHC, both 
in hard and soft processes, suggests that a new, strongly interacting state-of-matter has 
been created, the quark-gluon plasma (QGP). One of the most striking observables is the 
suppression of high-pj_ particles ("jet-quenching") in A-A compared to a simple scaling of 
p-p collisions [1-4], which together with dijet momentum imbalance [5, 6] finds a natural 
explanation in parton energy loss in the QGP [7]. Another observable in heavy- ion collisions 
which received much attention is heavy-quarkonium suppression. Such a suppression is 
expected from the Debye screening of the in-medium heavy-quark potential, and was thus 
originally proposed as a potential signal of QGP formation (and a direct probe of the 
plasma temperature) [8] . However it was later realized that several other effects can modify 
quarkonium yields in A-A collisions (see [9] for a review), some of those effects playing a 
role even in absence of any hot medium. 

In order to quantify the properties of the QGP created in heavy-ion collisions - a 
main goal of the forthcoming measurements at RHIC-II and LHC - a solid understanding 



of the nuclear modification of particle spectra in cold nuclear matter is thus required. A 
baseline for the study is provided by p-A (or d-A) collisions where a significant suppression 
is already reported for some particle species. In particular, light hadron [10, 11] and 
J/V> [12, 13] production in p-A collisions at forward rapidity is significantly lower than 
expected from the naive scaling of p-p spectra. 

At the same time it is astonishing that no consensus on the cold nuclear matter effects 
responsible for 3/i/j suppression has been achieved yet. In addition to the modifications of 
nuclear parton distribution functions, various mechanisms have been proposed to explain 
J/-0 suppression in p-A collisions. In the nucleus rest frame, a high-energy ,1/ip is formed 
long after the nucleus thus what actually propagates through the nucleus is the parent cc 
pair. Some approaches attribute 3/ip suppression to an effective absorption cross section 
Cabs of the cc pair (see Refs. [14] and [15] for recent works). Other models attribute J/ip 
suppression to the increase of the cc pair invariant mass by the multiple soft rescatterings 
through the nucleus, leading to a reduction of the overlap with the J/tp wave function [16] 
(see also [17]). In the approach used by [18] and that of our group [19], the dominant role is 
played by parton radiative energy loss, i.e., gluon radiation induced by multiple scattering 
of fast partons (or color octet cc pair) travelling through the nucleus. 

The present study, together with our previous works [19-21], supports parton energy 
loss induced by pj_-broadening as the main effect in J/ip suppression. At this point we 
should stress that the physical content of our approach and that of [18], both based on 
parton energy loss, are actually quite different. In [18], it is assumed that the interference 
between gluon emissions off the incoming and outgoing parton participating to the hard 
production process can be neglected. Under this assumption the induced energy loss of an 
incoming gluon in J/ip production is parametrically the same (up to color factors) as that 
of the incoming quark in the case of Drell-Yan pair production [22]. 

It was however argued in [21] that when the outgoing parton (or compact cc pair in 
the case of J/ip production) is produced at small angle (i.e., large energy E at limited 
p±) in the target nucleus rest frame, the medium- induced gluon spectrum is dominated 
by the interference between initial and final state radiation of gluons with large formation 
times. Such radiation is expected in low-p± J/i/j production in p-A collisions, where an 
incoming gluon is scattered at small angle into a compact color octet cc pair. The associated 
energy loss is proportional to the J/ip energy E and the role of this effect is thus expected 
to increase with increasing rapidity. This observation is at the basis of the energy loss 
scenario used in [19] to describe i/ip suppression as a function of y (or x F ) in a wide 
collision energy range for minimum bias p-A and d-A collisions. In the present study we 
generalize this approach to address the p±_ and centrality dependence of nuclear suppression 
at fixed-target (E866), RHIC and LHC energies. 

The paper is organized as follows. In Section 2 we generalize the calculation scheme 
of [19] to address the p± and centrality dependence of J/ip nuclear modification factors. 
Section 3 is devoted to the comparison of the model to E866 and latest PHENIX data and 
predictions in p-Pb collisions at the LHC are given in Section 4. Results are summarized 
in Section 5. 



2 Model 

2.1 Shift in E and p x 

Heavy-quarkonium 1 nuclear suppression in minimum bias p-A collisions as compared to 
p-p collisions can be represented in terms of the ratio 

z?V> i \ l dcr pA / dcr^p 

< (y '^ ) = i^dWdid^' {2A) 

where y and p± are the quarkonium rapidity and transverse momentum in the cm. frame 
of an elementary p-N collision (of energy y/s). By convention forward (positive) rapidities 
correspond to the proton fragmentation region. 

In the present study we generalize the model of Ref. [19] by expressing the quarkonium 
double differential cross section in p-A collisions in terms of that in p-p collisions, where 
a shift in the quarkonium energy E (defined in the nucleus rest frame) as in Ref. [19] but 
also a shift in p± account, respectively, for the energy loss e and transverse momentum 
broadening Ap± of the octet QQ pair propagating through the nucleus. The "double shift" 
in E and p± relates the p-A and p-p double differential cross sections da/dEd 2 p± as 2 

1 dfT pA /■/"», „, d<ri 



A dE d 2 pj 



LL v ^ e) A {e+ ^- a ^- (2 - 2) 



The quantity V(e, E) is the energy loss probability distribution or quenching weight (see 
Section 2.2) associated to the medium-induced radiation spectrum in a target nucleus A 
as compared to a lighter (e.g., proton) target nucleus. The integral over e is bounded by 
£max = min(-Ep — E, E), where E p ~ s/(2m p ) is the projectile proton energy in the nucleus 
rest frame. (We work in the limit y/s 3> m p , with m p the proton mass.) For the time being 
we consider minimum bias p-A collisions, hence the normalization factor 1/A in the l.h.s. 
of (2.2), but the model will be generalized in Section 2.5 to p-A and d-A collisions in a 
given centrality class. We will assume that Ap± is uniformly distributed in the azimuthal 
angle ip, and has a modulus defined by (2.10). 

The relation between the p-A and p-p differential cross sections in y and p±, 
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can be simply obtained from (2.2) by using 

In I 

E p M ± 



y(£,P±) = hi(##), (2.4) 



1 denoted by "ip" i n the rest of the paper. 

2 The dependence of Ap± on the azimuthal angle ip, Ap± = Ap±(ip), will be implicit in the following, 



and we will use the notations: 



— ; / = / de ; / = / d 2 p± etc. 



where M± = (M 2 + p\) 2 , with M the mass of the QQ pair. In Eq. (2.3), the variable 
E is given by E = E(y,p±) = E p e v M^/t/s, and the p-p cross section is evaluated at 
y = y{E + e,p± - Ap ± ). 

For the double differential p-p cross section da pp /dyd 2 pj_, we will use a parametrization 
consistent with the p-p data, see Section 2.3. In our model, the nuclear modification factor 
(2.1) obtained from (2.3) is thus fully determined by the quenching weight V(e, E) and the 
transverse momentum broadening Ap±. 

2.2 Quenching weight and Ap± 

In Ref. [19] the appropriate medium-induced gluon radiation spectrum dl/du in quarko- 
nium production in p-A collisions (as compared to p-B collisions) was derived, as well as 
the associated quenching weight V(e,E), 

dl ( f°° d/1 

V(e,E) = -ex P \- du— , (2.5) 
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where £j_a represents the accumulated transverse momentum transfer due to soft rescatter- 
ings in the target nucleus A, and Ab = max(A QCD ,£j_B). The latter dependence on A QCD 
arises from vetoing induced gluon radiation with k± < A QCD [19]. Note that V(e,E) is 
determined analytically in terms of the dilogarithm function Li2(x), 

V(s,E) = -exp \—Z-± Li 2 --L -Li 2 (-^-^) . (2.7) 
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(2.8) 

where La is the effective path length 3 across the target nucleus A and <?a the transport 
coefficient [19] 

. . /10- 2 \°- 3 . , 1 M± _„ . . 

QA = q{XA) = go ; ^A = min(3;oA,Z2) ; x 0A = „ — r~ ; x 2 = ^^e y . (2.9) 

\ x A J 2m p L A V s 

The parameter % = 0.075 GeV /fm was extracted in Ref. [19] from a fit to the E866 
data for Rj/} Be - We will use this value in the present study, which thus contains no free 
parameter. 4 



3 For the present study we will use Lb c = 3.24 fm, Lf c = 6.62 fm, Lw = 9.35 fm, La u = 10.23 fm and 
ipb = 10.11 fm for minimum bias collisions, as well as L p — 1.5 fm [19]. The dependence of Lau (resp. 
Lpb) on the centrality class in d-Au collisions at RHIC (resp. p-Pb collisions at LHC) is described in 
Section 2.5, see Table 2, and in Appendix B. 

4 Let us also note that in the fitting procedure of Ref. [19], small values of qo and thus of qL were 
explored, for which the dependence of (2.6) on A QCD is relevant. However, it turns out that the extracted 
go = 0.075 GeV 2 /fm is large enough to satisfy \/q L > 0.25 GeV for all values of La (including L p ) and xi 
considered in the present study, as can be easily checked from (2.9) using L p — 1.5 fm. The present study 
is thus independent of the value of A QCD provided A QCD < 0.25 GeV. 



Finally, the transverse momentum broadening in p-A with respect to p-B collisions is 
simply given by 

(Ap ± ) 2 = £ 2 ±A - £ 2 ±B = q A L A - q B L B . (2.10) 

In our study we neglect the fluctuations of the broadening i\ around the average value 
qL. We have checked that assuming a distribution in l± of the Gaussian type, P{£±) oc 
exp (—i^/qL), only slightly modifies (by at most 10%) the predictions for Rz A (y,p±) pre- 
sented in Sections 3 and 4, without changing the overall shape of i? p A- 

2.3 Parametrization of the p — p cross section 

Similarly to Ref. [19], we use for the double differential p-p cross section da pp /dyd 2 p± a 
simple parametrization consistent with the available p-p data, rather than relying on some 
model-dependent quarkonium production mechanism in hadronic collisions. 

The double differential cross section of prompt J/i/j and T production can be conve- 
niently parametrized as 

d °t - Mx (PL\ m X ( l _^ COsh y) n ^ MxKp±)xH y, p±) . (2.11) 



dyd 2 p± \pI+p\ 

At the LHC {yfs = 7 TeV) the values of the free parameters jV, po, m and n are obtained 
from a global fit of ALICE [23], ATLAS [24] and LHCb [25] data on prompt J/t/> production 5 
and from a fit of LHCb data [26] on T production. They are summarized in Table 1 together 
with the corresponding xV n df values. 6 



Quarkonium 


v^ (GeV) 


Po (GeV) 


n 


m 


X 2 /ndf 


T 


7000 
7000 


4.2 
6.6 


19.2 
13.8 


3.5 

2.8 


54/139 

48/72 


J/V> 


200 
38.7 


3.3 

3.1 


(8.3) 
(4.5) 


4.3 
5.3 


27/37 
45/19 



Table 1. Values of the fit parameters obtained from a fit to LHC (J/ip and T), PHENIX, and 
E789 p-p data. Values given in parenthesis are fixed in the fitting procedure (see text for details). 

At RHIC the amount and precision of data are not sufficient to fix precisely the fit 
parameters. The fit to the double differential data measured by PHENIX [27] is therefore 
performed by fixing the value of n = 8.3 obtained from the fit of the single differential cross 
section da pp /dy, see Ref. [19]. 

Finally, a fit to the J/tp E789 (y/s = 38.7 GeV) data [28] is also performed with the 
aim to compare the model predictions to the E866 data [12] at the same center-of-mass 
energy. 7 Like at RHIC, the number of fit parameters is reduced by fixing the value of 
n = 4.5 obtained from fitting da pp / dx F data at this energy [19]. 



5 The ALICE data are only given for inclusive 3/t/j production, i.e. including 3/tp from B-decays. 
6 The value of A/" is irrelevant as we consider only cross section ratios and is therefore not given. 
7 Note however that the E789 data have been taken at x F = only. Therefore this parametrization may 
not be valid at large x F for which the model is compared to the E866 measurements, see Section 3.1. 
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Figure 1. Comparison between heavy-quarkonium (J/i/>, T) production data in p-p collisions and 
the parametrization (2.11) (solid red line). Data arc taken from ATLAS [24], LHCb [25, 26], and 
PHENIX [27]. 

The agreement between some of the prompt J/i/j and T measurements at RHIC and 
LHC and the parametrization (2.11) is shown in Fig. 1. Finally, let us mention that the 
functional form of the parametrization (2.11) is fully consistent with that for the single 
differential rate in x F used in [19], as we briefly show in Appendix A. 

2.4 A useful approximation for -Rp A (y,p_i_) 
Using (2.3) and (2.11) the attenuation factor (2.1) reads 



#p A (?/>P-U 



V(s,E) 



ifi Je 

V(e,E) 
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E + e 
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Since V(e, E) is peaked at small values of e, we neglect e in the latter (p integral. In this 
approximation, the ip and e integrals factorize, 



loss/ 
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The factor R^™ ad (y,p±) describes nuclear modification due to transverse momentum 
broadening only, as can be seen by setting V(e, E) = 5(e) in (2.12). The factor R l °^(y,p±) 
describes the effect of energy loss only, obtained by setting Ap± = in (2.12). In the 
following we will use the factorized expression (2.13), which turns out to be a very accurate 
approximation to (2.1) in all the practical applications of Sections 3 and 4. 

Finally, let us mention that the p^-inclusive nuclear suppression factor -Rp A (y) studied 
in Ref. [19] can be simply recovered from (2.12), along the same lines as in Appendix A 
where the parametrization of da pp /dy is obtained from d<7 pp /dyd 2 pj_. Integrating both the 
numerator and denominator of (2.1) (and thus of (2.12)) over p±, the function v(y,pi_) can 
be replaced by its value at a typical p± determined by the width of fJ>(p±) (see (A. 6)), and 
the pj_-integral of /x(p±) cancels between the numerator and denominator. As a result, the 
p_i_-inclusive suppression factor reads R p a(.u) — ^ p ° A s (?/,p_i_), which corresponds exactly to 
the quantity studied in Ref. [19]. 

2.5 Centrality dependence 

In the preceding sections we addressed the case of minimum bias p-A collisions. The 
model is generalized to the case of p-A collisions at a given centrality (or at a given impact 
parameter b) by using the effective length La corresponding to that centrality. Since La 
fully determines the essential quantities £±a and Ap± of the model (see (2.8), (2.9) and 
(2.10)), this is the only modification required. However, the impact parameter (as well as 
the number of participants) is not a direct experimental observable, and the consistent way 
for making a theoretical estimate of La is to follow the experimental procedure as closely 
as possible. 

Both at RHIC [29] and the LHC [30] a centrality selection is done by triggering on 
event multiplicity in forward detectors. This multiplicity is strongly correlated with the 
number of participating nucleons from the target nucleus. Thus, multiplicity cuts impose a 
restriction on the number of participants in a given event. The values of the forward event 
multiplicities which separate centrality classes in the experiment are chosen to attribute a 
certain fraction of the total inelastic cross section to each centrality class. The common 
choice is making 20% slices from the most central (largest multiplicity) to the peripheral 
(lowest multiplicity) events. The exception is the most peripheral class which is taken at 
60-88% of the total inelastic cross section at RHIC (class D) and 60-100% at the LHC 
(class 4). 

To compute the average path length La for each centrality class, we employ the fol- 
lowing procedure. 

First, we define centrality classes in terms of the number of participants N p . Within 
a Glauber description (see Appendix B.l), we define a centrality class by the threshold 
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Armin. Armax 
Jv p 5 Jv p 


P(class) 
P(JV>1) 


(N c ) 


Lau 


class 


7\rmin . lumax 


P(class) 
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(No) 


Lpb 


A 


11; 197 


0.28 


15.9 


12.87 


1 


12; 208 


0.246 


14.8 


13.46 


B 


8; 12 


0.24 


10.9 


9.62 


2 


9; 12 


0.215 


10.5 


9.55 


C 


5; 8 


0.23 


7.0 


7.17 


3 


5; 8 


0.215 


6.5 


6.29 


D 


2; 4 


0.29 


3.6 


3.84 


4 


1;5 


0.428 


2.4 


3.39 



Table 2. Average number of binary collisions (N c ), and average path length La (in fm) for different 
centrality classes (A-D at RHIC and 1-4 at LHC), as obtained in a Glauber calculation. See text 
and Appendix B for details. 

values iV™ m and N™ 3 * of N p which saturate approximately the same fraction of the total 
interaction probability P(N P > 1) = 1 as the fraction of the total inelastic cross section 
attributed to the centrality classes in the experimental selection procedure. We note how- 
ever that centrality in the experiment is defined in terms of event multiplicity rather than 
N p . In a given multiplicity class one may have events where N p is slightly above or below 
the thresholds defined according to a sharp cut on the fraction of the total interaction 
probability as a function of N p . In order to account for this possibility we widen the N p 
interval attributed to each class and use the new threshold values in further estimates. 
Those values are given in Table 2. As a consistency check of our class selection method, we 
calculate the average number of binary collisions (N c ) for the different centrality classes. 
These numbers coincide with the results of the Glauber Monte-Carlo supplemented with 
the p-p RHIC data used by the PHENIX collaboration (see Table 1 in Ref. [13]). 

Second, within each of the centrality classes defined in such a way, we determine, also in 
the Glauber model, the average number of target nucleons participating in the rescattering 
of the fast color octet QQ pair, see Appendix B.2. The average path length La for a 
given centrality class directly follows from this number. Table 2 displays the values of L\ u 
(RHIC) and Lpb (LHC) to be used in our model, for each centrality class. 



3 Comparison to E866 and RHIC data 

The J/ip nuclear production ratio is computed as a function of p± for various values of y 
or x F = (2Mj_/i/i) sinhy. The only parameter of the model, the transport coefficient %, 
has been fixed to qo = 0.075 GeV 2 /fm in [19] from the x F dependence of J/ip suppression 
measured in p-W collisions by E866 [12]. In the numerical calculations we use A QCD = 
0.25 GeV (however see footnote 4), a s = 0.5, and M = 3GeV (M = 9GeV) for the mass 
of the cc (bb) pair. 

3.1 E866 

The E866 collaboration has measured the J/^ suppression in p-Fe and p-W collisions 
(with respect to p-Be) at y/s = 38.7 GeV as a function of the transverse momentum 
for three domains in x F [12]. At small x F , J/tp production can be affected by nuclear 
absorption since the typical J /tp hadronization time becomes comparable to (or less than) 
the size of the nuclear targets [19]. The model is therefore compared to the E866 data 



in the intermediate-x F (0.2 < x F < 0.6, (x F ) = 0.308) and large-x F (0.3 < x F < 0.93, 
(x F ) = 0.48) domains. 



E866 Vs = 38.7 GeV — E. loss + broad. - - Broad 
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P T (GeV) 
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Figure 2. Model predictions (solid red curves) for the J/i/" nuclear suppression factor compared to 
the E866 data for i?Fo/Bo(pj-) (left) and Rw/Bo(p±) (right), in the intermediate- x F ((x F ) = 0.308, 
top) and large- a;^ ((x F ) = 0.48, bottom) ranges. The dashed lines indicate the effect of momentum 
broadening only, Rl^ ad (y,p ± ), Eq. (2.14). 

In Fig. 2 are shown as solid lines the Fe/Be (left) and W/Be (right) nuclear production 
ratios at intermediate (top) and large (bottom) x F . The ratio R p a(p±) increases withpj_ in 
almost the whole p± range, with more pronounced suppression (at pj_ = GeV, where R p a 
is the smallest) in W targets and at large x F . The pj_ dependence essentially arises from 
that of Rp™ ad (p±), shown as dashed lines. The sole energy loss effect, R l " A s {y,p±), proves 
rather flat in this pj_-domain, but is essential to fix the magnitude of -R p a, leading to a 
remarkable agreement between the data and the model predictions in the < p± < 2 GeV 
range. 

At large x F , the data overshoot the theoretical expectations above p± > 3 GeV, for 
which the model predictions flatten out. Several reasons might explain this disagreement. 
First of all, the parametrization of the p-p cross section (using E789 data, see Section 2.3) 
has been performed at x F = only; it could therefore well be that this fit, used to compute 
i? p A, is no longer appropriate to describe the p-p production cross section at both large 
x F ~ 0.4 and p± > 3 GeV. Also note that the theoretical calculations have been performed 
at x F = {x F ) while the data are averaged on a rather large x F -bin. When p±_ gets larger, 
the kinematical correlation between p± and x F would therefore tend to decrease the typical 



x F , for which the model would predict slightly larger R p a ratios. 8 For those reasons it is 
difficult to draw any firm conclusion from the comparison between the model and the E866 
data at large x F and p± . 

Putting aside the latter (p±,x F ) region, it is remarkable that the model reproduces 
quantitatively the p± dependence of J/ip suppression in various nuclei and at different x F 
values. The fact that the same quantity, qL, determines the strength of medium- induced 
gluon radiation (and therefore energy loss) necessary to explain the x F dependence of R p a 
(see Ref. [19]) and the amount of momentum broadening required in the present study to re- 
produce the shape of R p \ as a function of p± , strongly supports parton energy loss induced 
by momentum broadening as the dominant effect in quarkonium nuclear suppression. 

3.2 RHIC 

Let us now move to RHIC energy, where the p± dependence of J/ip suppression in d-Au 
collisions has been reported recently by the PHENIX collaboration [13]. 

In Fig. 3 the model predictions are compared to the PHENIX data measured in mini- 
mum bias d-Au collisions, at backward (—2.2 < y < —1.2, left), central (\y\ < 0.35, middle) 
and forward (1.2 < y < 2.2, right) rapidities. 9 The model reproduces the trend seen in 
data, namely an increase with p± up to p± ~ 4-5 GeV. Around those values, some nuclear 
enhancement, R p a(p±) > 1, is visible at backward and central rapidities, but would need 
more precise data to be confirmed. At forward rapidity, the suppression due to energy loss 
is too strong to observe such an enhancement, both in the data and in the model. As for 
the results at E866 energy discussed earlier, the p± shape of R p a(p±) is essentially driven 
by the effect of momentum broadening, Eq. (2.14), shown as dashed lines. 

On top of minimum bias collisions, J/ifi suppression has also been measured in four 
centrality classes of d-Au collisions (0-20%, 20-40%, 40-60%, 60-88%). The data are 
shown in Fig. 4 for the three rapidity bins and four centrality classes and compared to 
the model. Details on the medium length L corresponding to the various centrality classes 
of d-Au collisions and used in the theoretical calculation can be found in Section 2.5 and 
Appendix B. As for minimum bias collisions, the model predictions prove in very good 
agreement with data. In particular, the centrality dependence is well reproduced by the 
model, with rather pronounced effects in the 0-20% most central collisions, and almost 
negligible effects, R p a — 1 in the whole pj_-range for the 60-88% most peripheral collisions. 

Therefore it appears that the present model, based on parton energy loss and momen- 
tum broadening, offers a better agreement on the p±_ and centrality dependence of R p a than 
models based purely on modification of parton densities [15] and nuclear absorption [14] 
which tend to predict a flatter dependence of R p a(p±)- 



We checked that averaging properly over the whole x F -range considered experimentally increases by 
~ 5% the 7?pA ratios at large p±, therefore slightly reducing the discrepancy with experimental data. 

Calculations are performed at the fixed values of y — —1.7, 0, +1.7 respectively. We checked that 
similar results are obtained when averaging over the experimental y range. Note also that for y — —1.7, 
the J/ip hadronization time is comparable to the size of the gold nucleus, and nuclear absorption may play 
a role [19]. 
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Figure 3. Model predictions for the J jt\> nuclear suppression factor R p a(p±) in minimum bias d-Au 
collisions at RHIC, at backward (left), central (middle) and forward (right) rapidities (solid curves). 
The dashed lines indicate the effect of momentum broadening only, R^ a ' d {y,p±), Eq. (2.14). 
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Figure 4. Same as Fig. 3 in the four centrality classes (from top left to bottom right in each panel) 
at backward, central and forward rapidities (upper left, upper right, bottom). 
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4 Predictions for p— Pb collisions at the LHC 

A run of p-Pb collisions at y/s = 5 TeV has taken place at LHC early 2013 with an 
integrated luminosity of C ~ 30 nb _1 collected by ALICE, ATLAS and CMS experiments. 
This will allow for precise measurements of J/?/> and T production in p-A collisions at an 
unprecedented energy, on a wide range of rapidities (\y\ < 5) and transverse momenta, and 
consequently clarify the role of cold nuclear matter effects at high energy. In this section, 
we therefore provide predictions for J/ip and T nuclear production ratios R p a as a function 
of p±, for different rapidities (y = —3.7, 0, 2.8) 10 and centrality classes (labelled 1. . . 4 from 
central to peripheral collisions, see Section 2.5 for details). 

In Fig. 5 we show the J/?/> nuclear production ratio for minimum bias p-Pb collisions, 
from backward (left) to central (middle) and forward (right) rapidities. At all rapidities 
a depletion of J/ip production is expected at low p±, say p± < 3 GeV. At larger p± a 
"Cronin peak" might only be visible at large rapidity, y = 2.8 and above, in minimum bias 
collisions. Quite generally the effects of both energy loss and momentum broadening are 
expected to become more pronounced at larger rapidities. As shown in [19], energy loss 
effects prove stronger at large positive rapidity because of the energy dependence of the 
average energy loss, AE oc E, associated to the medium-induced spectrum (2.6). 






p T (GeV) 



p T (GeV) 



Pt (GeV) 



Figure 5. Prediction for the J/^ nuclear suppression factor R p a(p±) in minimum bias p-Pb 
collisions at LHC, for backward, central and forward rapidities. 

The model predictions are also provided in the four centrality classes of p-Pb collisions 
in Fig. 6. As expected the deviations of R p \ from unity are largest in the most central 
collisions, while in the most peripheral p-Pb collisions (centrality class 4), R p a(p±) — 1 
at all p± > 2-3 GeV. The most spectacular effects can be seen in the central collisions 
(class 1) and at forward rapidity, where R p a — 0.25 at p± = GeV and R p a — 1-3 at 
p± = 6 GeV. 



10 Here the rapidity y is defined in the center-of-mass frame of the proton-nucleon collision, related to the 
rapidity in the laboratory frame ?/i a b as y — j/i a b — 0.465 in p-Pb collisions and y = j/i a b + 0.465 in Pb-p 
collisions. The value y — 2.8 (respectively y — —3.7) is chosen to correspond to the median rapidity of the 
ALICE spectrometer acceptance, namely 2.5 < yiab < 4, in p-Pb (respectively Pb-p) collisions. 
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Figure 6. Same as Fig. 5 in the four centrality classes (from top left to bottom right in each panel) 
at backward, central and forward rapidities (upper left, upper right, bottom). 



Predictions are performed as well in the T channel. The expected Bs A at the LHC 
is shown in Fig. 7 for the most central p-Pb collisions (centrality class 1). Because of 
the mass dependence of the energy loss, AE oc M7 , the suppression due to energy loss, 
-^pA 3 ' ^ s m ilder than for J/tp. Moreover, although the amount of momentum broadening 
experienced by J/tp and T states is expected to be similar (neglecting the x dependence 
of the transport coefficient), the height of the Cronin peak is much less pronounced for T 
than for J/ip (compare e.g. Fig. 6 bottom and Fig. 7 right). 11 This is due to the flatter 
T pj_-spectrum as compared to that of J/ip production, see Fig. 1 and the values of the 
parameters po an d m in Table 1. As can be seen in Fig. 7, T suppression, Bl A < 1, 



is 



predicted in the range < p± < 6 GeV at mid-rapidity. The suppression extends to larger 
p± than for J/tp due to the larger value of the po parameter in the p-p cross section. The 
suppression is maximal at p± = GeV, where Rl A ~ 0.85 (resp. 0.65) at y = (resp. 
2/ = 2.8). 



xl As a matter of fact, the slight enhancement arising from R p ™' 
R Y pA < 1, making R pA smaller than one at all p±. 



is compensated by energy loss effects, 
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Figure 7. Predictions for the T nuclear suppression factor R p a(p±) in central p-Pb collisions 
(centrality class 1) at LHC, for central (left) and forward (right) rapidities. 



5 Conclusion 

Following our earlier work [19], we studied the effects of parton pj_-broadening and energy 
loss in cold nuclear matter on the p± dependence of J ftp and T suppression in p-A collisions. 
We found that the momentum broadening is responsible for the fast variation of J/ip 
suppression with p±, while medium- induced energy loss essentially affects the magnitude 
of R pA . 

Using the transport coefficient % = 0.075 GeV 2 /fm fixed in [19], the model predic- 
tions prove in very good agreement with recent PHENIX data [13] in minimum bias and 
centrality-dependent d-Au collisions at \fs = 200 GeV. Our results are also successfully 
compared to earlier results from the E866 collaboration [12]. Finally, predictions for 3/tp 
and T suppression in p-Pb collisions (minimum bias and in four centrality classes) at the 
LHC (y/s = 5 TeV) are provided. 

The good description, within a consistent framework, of both the rapidity and trans- 
verse momentum dependence of Rjp[ from fixed-target experiments to RHIC is a hint 
that parton energy loss induced by momentum broadening might be the dominant effect 
responsible for J/i/> suppression in p-A collisions. 
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A Fits of p— p single and double differential rates 

Here we compare the parametrization (2.11) for the double differential p-p cross section 
used in the present study and that for the single differential rate d<T pp /dx F used in [19], 

< K (1^, (A.1) 

dx F x 

where the variables x F and x' are defined by 

E E p M? 2M X , . . . 

x * = e--^ — = ^7T smhy ' (A ' 2) 

. E E v M? 2M ± , , , 

x' = — + -^ ^ = —L coshy . A.3 

E p E s y/s 

Eq. (A.l) translates in rapidity to (use dx F /dy = x') 

^ oc (1 - x'T = (l - ^ coshy)" = i/(y, PjL ) . (A.4) 

In Ref. [19] no information on ^-distributions was used, and p± was replaced by some 
typical value, assumed to be p± = 1 GeV. 

In the present study, the single differential cross section da pp /dy can be obtained by 
integrating (2.11), 



data f da 



pp _ / ^"pp 



4> 



dy Jp, dyd 2 p ± 



N / v(p±)v(y,p±) -N 



m(pj 



;-'_ 



v(y,p±), (A.5) 



where the latter approximation arises from p(p±) decreasing much faster than v(y,p±) 
with p±. This can be checked for the values of the parameters po, m, n (see Table 1) and 
within the intervals in p± and y considered in our study (see Sections 3 and 4) . Comparing 
(A.4) and (A.5) we see that the parametrization (2.11) is consistent with that for the single 
differential rate (A.4) (or equivalently (A.l)) used in [19]. 

Let us remark that in (A.5) the typical p± may be defined as 

^ dp 2 ± fi(p ± ) = -] dp\ /i(p ± ) => pi = (^2 m - 1 - lj pi (m > 1) . (A.6) 

Using the values of po and m given in Table 1 for J/V', we find that p±_ varies in the 
range p±{po,Tn) = 1.3 — 2.4 GeV, somewhat above the ad hoc value p± = 1 GeV used 
in [19] . However, taking p± ~ 2 GeV instead of 1 GeV affects only slightly the value of M± 
(M± ~ 3.6 GeV instead of 3.2 GeV), with no sizeable effect on the predictions presented in 
Ref. [19]. 
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B Effective path length vs. centrality class 

In this Appendix we explain how the numbers of Table 2 in Section 2.5 have been obtained. 

B.l Number of participants and binary collisions for a given centrality class 

p— A collisions (LHC case) 

In the Glauber model, at a given impact parameter b the number of participating nucleons 
N p in the target nucleus follows a binomial distribution, 

Jdb(A)l^(b)l^i-mW]'-* 

/db{i-[i- PpA (b)]- 4 } 

Here p p \(h) = u^ ?A.(b) is the probability for an inelastic collision between the projectile 
proton and a nucleon of the target nucleus, and <r^ N is the inelastic p-p cross section 
which we take to be 70 mb at LHC energies. The target nucleus optical thickness Tj^(h) is 
normalized as JdbT A (b) = 1. The denominator in (B.l) is the p-A total inelastic cross 
section af and ensures the correct normalization P(N P > 0) = 1 for the total interaction 
probability. 

The number of binary collisions coincides with the number of participating nucleons 
of the target. 

d— A collisions (RHIC case) 

The generalization of Eq. (B.l) to the d-A case is straightforward, 

/db/drP d (r)(^)[p (U (b,r)]^[l-J»dA(b,r)] A - Ar '' 

p(jV ) = - - p t (B.2) 

Jdb/drP d (r) {l - [1 -p dA (b,r)] A } 

and is obtained by introducing the distribution -Pd( r ) f° r the p-n transverse separation 
r in deuterium, and replacing p pA (b) by the interaction probability p dA (b,r) of a target 
nucleon with either nucleon of the deuterium projectile, given by 

p dA (b, r) = o£N T A (b + §) + a™ T A (b - §) - P2 (b, r). (B.3) 

Here 

p 2 (b, r) = J dsT A (s)p(b + r/2 - s)p(b - r/2 - s) ~ T A (b) f dsp(r + s)p(s) (B.4) 

is the collision probability of a target nucleon with both nucleons of the deuterium projec- 
tile. The denominator of (B.2) is the d-A total inelastic cross section af n . 

The probability (B.4) depends on the N-N inelastic collision probability as a function 
of the impact parameter p(s). For that profile at cm. energy y/s = 200 GeV we take 
the Regge-inspired parametrization p(h) = 1 — exp(— 2Ne~ h ' a ) with a = 1.05 fm 2 and 
N = 1.1, giving the total p-p inelastic cross section o"^ N = Jdsp(s) = 4.2 fm = 42 mb. 
The distribution -P d (r) is evaluated by assuming a Hulthen form for the deuterium wave 
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function. The thickness functions of all target nuclei considered in our study (Be, Fe, W, 
Au, Pb) are extracted from low-energy electron-proton scattering experiments [31]. 

Summing up the probabilities (B.l) and (B.2) we find the threshold values N™ m and 
jymax w ] a j c j 1 saturate approximately 20% of the total interaction probability P(N p > 1) = 
1, as described in Section 2.5. 

In the d-A case, the number of binary collisions differs from N p by the number of 
target nucleons which undergo collisions with both nucleons of the deuterium projectile: 

N 2 

E N( A ) j db / drP d (r)p 2 (b, r) [p dA (b, r)] N ~ x [1 - p dA (b, v)] A ~ N 

(AN)\ Ne[NuN2] = 2—^- . 

E (£) fdhf drP d (r) [p d A(b, r)]" [1 - PdA (h, r)] A ~ N 
N=Ni 

(B.5) 

B.2 Average number of collisions in the events with a hard process 

p A (LHC case) 

In the case of triggering on J/ip production (or any other hard process with a small 
cross section) the average number of collisions for the participating nucleon gets mod- 
ified. The hard production process can occur in each of the inelastic collisions of the 
projectile nucleon with a probability aj/^/a™ , 12 so that the hard process probability 
is P(J/i/j\N p ) = Npaj/^p/a^ 1 and the joint probability for the hard production and N p 
participating nucleons in the target is 

P(N P , J/VO = ^^ fdb( A ) [p pA (b)]^ [1 - PpA (h)] A - N * . (B.6) 

The normalized probability distribution for the number of collisions of the projectile 
proton in the events tagged by both centrality and J/ip production is 

3&iVp /db(£) [PpA(b)]^ [1 - PpA (b)] A ~ N v 

P(N p \J/4>, [N U N 3 ]) = -Jj 1 (B.7) 

E ^Nf db( A ) [ PpA (h)] N [1 - PpA (b)] A - N 

N=N 1 in 

and is independent of the hard process cross section. The corresponding average of N p is 

E N(N - 1) / db(jj) [ PpA (h)] N [1 - PpA (b)] A - N 

(N P )\ J/nNl , m] = 1 + ^^ • (B.8) 

E N J dh{ A ) [p pA (b)] N [1 - PpA (h)] A - N 

N=Ni 

The interpretation of (B.8) is straightforward. Unity stands for the target nucleon which 
participated to the hard process. The second term corresponds to the target nucleons 



12 We assume that the hard process cross section is small, so that the probability of two hard processes 
in the same p-A collision can be neglected. 
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which also undergo an inelastic collision with the projectile, and may thus contribute to 
the transverse momentum broadening of the cc pair, with the probability Abroad/ °to- 

The effective path length in the target nucleus for the centrality class [A^i,AT 2 ] thus 
reads 

N 2 

E N(N - l)/db(J) [p pA (b)f [1 - PpA (b)] A ^ 
La = L p+ N=N ^ , (B.9) 

<%o E N f dh(%) [ PpA (h)] N [1 - PpA (b)] A - N 

where L p is the corresponding length in a proton target. For the minimum bias case doing 
the summations is trivial and one recovers the expression used in Ref. [19]. In the numerical 
applications we take L p = 1.5 fm and po = 0.17 fm~ 3 . 

d A (RHIC case) 

Compared to the p-A case some complications arise because of the deuterium projectile. 
What is relevant for the broadening is the number of collisions Nt suffered by the nucleon 
of the deuterium participating to the hard process. This number is not equal to the number 
of participants in a given event. In the binomial expansion of the first multiplier in the 
numerator of (B.2) 

[pdA(b,r)]^ = ($) \p pA (b + r/2)] Nt \p pA (b-r/2)-p 2 (b,r)] N *- N * , (B.10) 

each term corresponds to the probability of N t collisions of the tagged nucleon of the deu- 
terium projectile with overall N p participants in the target nucleus. The joint probability 
of the production process in d-A collision with given Nt and overall N p thus reads: 

2N 



P(N p ,N t ,J/i/>) = -^^ fdb [drP d (r)x 

"in ^in J J 



in in 



(4.) @ [P P A(b + r/2)]"« [p pA (b - r/2) - p 2 (b, r)] N '~ N * [1 - PdA (b, r)] A ~ N > . (B.ll) 

Correspondingly the number of inelastic rescatterings for the nucleon of the deuterium 
projectile tagged by the hard process is 

( JV *)lj/^,[JVi,JVa] = 1+ 

N 2 

E N(N - 1) f db f drP d (r)(£)b pA (b)] 2 [p dA (b,r)] N - 2 [1 - p dA (b, r)] A ~ N 
E N(N - 1) / db / drP d (r) ( A ) PpA (b) \p dA (b, r)]^ 1 [1 - PdA (b, r)] A ~ N 

N=Ni 

Interpretation of (B.12) goes along the same line as for (B.8). 
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The path length of the cc pair for a given centrality class is thus (after substituting 
p pA (b)=a™T A (b)): 

£ N(N - 1) / db / drP d (r) (£) (T A (b)) 2 [ PdA (b, r)]^ 2 [1 - p dA (b, r)] A "" 

L A = L P + 

Po £ iV(JV-l)/db/drP d (r)(5)r A (b)[p dA (b,r)] Ar - 1 [l-p dA (b,r)] A - JV 

JV=JVi 

(B.13) 

For the minimum bias case (N\ = 1, N2 = A) the summations in (B.13) are explicit and 

one again recovers the expression (3.18) of Ref. [19]. The path length for the different 

centrality classes is given in Table 2. 
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